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SECTION  I 

APERTURE  REPRESENTATION  BY  A PAIR  OF  DIPOLES 

A convenient  approximation  of  the  electronagnetic  field  in  the 
presence  of  a small  aperture  is  in  terms  of  a pair  of  dipoles,  as 
originated  by  Bethe  [1] . Assume  that  the  local  distribution  of  the 
total  electric  field  over  an  aperture  in  x-z  plane,  such 

as  in  Fig.l.(a),  has  been  determined  by  an  analytical  or  a numerical 
solution  of  the  boundary-value  problem.  Then,  the  electromagnetic 
field  in  the  "internal"  region  y>0  remains  unchanged  if  the  aperture 
is  closed  by  a metal  lid  on  top  of  which  there  is  a magnetic  surface 


current 


5g"*(x,z)  = t^(x,z)  X 


v^ere  a^.  is  the  unit  vector  in  y direction. 

The  coupling  to  the  T^l  wave  on  the  system  of  conductors  can 
be  then  computed  by  replacing  the  aperture  with  a conducting  lid  on 
top  of  which  there  are  two  dipoles,  and  c^^^  as  shown  in  Fig.  l.(b). 
Electric  dipole  moment  is  oriented  in  y direction; 


V 'jj 


and  the  magnetic  dipole  moment  of  interest  here  is  oriented  in  x direction; 


‘^mx  * 


z)dxdz. 
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The  magnetic  dipole  moment  may  also  have  the  z-canponent,  but 
the  present  report  is  devoted  solely  to  a coupling  of  quasi-TEM  waves 
guided  along  z direction  by  a system  of  conductors  as  indicated  in 
Fig.  l.(b).  Therefore,  the  z component  of  the  magnetic  dipole  is 
of  no  importance,  because  it  does  not  interact  with  these  quasi-TEM 
waves . 

In  this  report,  the  electric  dipole  moment  c^  (in  Ampere ‘meters) 
has  a meaning  of  the  moment  of  an  electric  current  element,  similarly 
as  in  references  [2] -[4].  This  current  moment  should  not  be  confused 
with  the  electric  charge  moment  p (in  Coulcmib  • meters)  such  as  used 
for  example  in  references  [5] -[7].  The  relationship  between  these 
two  moments,  for  exp(ja)t)  variation,  is  as  follows: 

• (4) 

Similarly,  the  magnetic  dipole  moment  c^^^  (in  Voltameters)  in 
this  report  denotes  the  moment  of  the  magnetic  current  element,  in 
the  saise  as  used  by  references  [8] -[10],  The  magnetic  charge  moment 
iff  (in  Ampere* square -meters),  such  as  used  for  instance  in  references 
[5] -[7]  is  related  to  as  follows: 

Cjjj  = jojimi.  (5) 


Figure  1. (b)  is  a first-order  equivalent  of  the  original 
configuration  from  Fig.l.(a),  Sometimes,  it  is  convenient  to  further 
change  Fig.  l.(b)  into  an  equivalent  configuration  in  Fig.  l.(c)  in 
which  there  is  a distribution  of  the  surface  magnetic  current  3*^  and 
of  the  surface  electric  current  over  the  x-y  plane.  Figures  l.(b) 
and  l.(c)  are  equivalent  if  the  surface  currents  become  delta  functions 


as  follows: 
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(6) 

(7) 


For  an  aperture  of  general  shape,  integrations  (2)  and  (3)  are 
to  be  performed  numerically.  For  several  characteristic  shapes  (circle, 
ellipse,  naiTOw  slit,  square,  etc.),  the  dipole  moments  have  been  com- 
puted or  determined  experimentally.  It  is  customary  to  express  the 
moments  in  terms  of  the  excitation  fields  2 , {?  and  in  terms  of  the 
polarizabilities  a^,  Consider  the  aperture  in  Fig.  2 vdiich  is  to 
be  replaced  by  the  dipole  moments  so  that  the  field  in  the  "internal” 
region  y > 0 is  maintained.  The  excitation  fields  (so-called  short 
circuit  fields)  produced  by  the  sources  located  in  the  internal  region 
will  be  denoted  and  When  the  excitation  cones  from  the 

side  y<0,  the  short  circuit  fields  will  be  denoted  by  and 

For  the  purpose  of  coupling  to  the  quasi-lEM  waves  in  the  internal  region, 
only  the  components  and  are  of  interest.  Instead  of  per- 
forming the  numerical  integrations  fron  (2)  and  (3),  the  dipole  moments 
of  small  apertures  of  characteristic  shapes  may  be  computed  as  follows: 

V ■ • t«) 

The  electric  polarizability  and  the  magnetic  polarizability  of  a 
circular  aperture  of  diameter  d are: 

8 


a-.r 


1 j3  1 ,3 

“e  ' rz  ' “m  " ^ • 


For  a square  aperture  of  side  I,  Cohn  [11],  [12]  has  measured  the 
following  polarizabilities: 

Og  = 0.1137£^  , = 0.2590JI^  . 


The  last  two  references  also  contain  the  measured  polarizabilities  for 
rectangular  and  other  shapes  of  apertures. 

In  Fig.  5,  a plane  wave  is  shown  incident  from  the  y<0  region. 
The  two  possible  polarizations  of  the  plane  wave  are  denoted  TE  and 
TT-l.  The  corresponding  dipole  mcments  are  computed  from  (8]  and  (9)  by 
substituting  the  following  excitation  fields; 


C ^ sin  a 


♦ 2Eo 

H ^ cos  e cos  a. 

xs  n 


For  this  excitation,  take  = 0 and  = 0- 

When  electromagnetic  fields  in  the  presence  of  a ground  plane 
are  considered,  sometimes  it  is  convenient  to  remove  the  plane  and 
replace  it  by  appropriate  images.  In  the  present  report  the  ground 
plane  has  not  been  removed  and  the  images  have  not  been  invoked. 


» 


Port  3 


TE  wave; 


I 


Fig.  3.  Orientation  of 

11 


I 


When  dipole  moments  and  have  been  determined,  it 
is  possible  to  canpute  the  amplitudes  of  the  outgoing  guided  waves 
created  by  these  dipole  manents.  In  Fig.  3,  the  scattering  ampli- 
tudes of  the  TBl  guided  waves  traveling  in  the  positive  and 
negative  z direction  are  denoted  by  a^  and  b^.  Using  the  appro- 
priate boundary  conditions  in  the  xy  plane,  the  following  scattering 
amplitudes  a^  and  b^  are  obtained  [13] : 

S • I VyTEM'^0-'”' ' 

*’s  * J 

The  scattering  amplitudes  are  complex  numbers,  normalized  in  such 
a way  that  lag|^/2  and  Ib^l^/Z  are  powers  of  the  outgoing  waves 
traveling  in  positive  and  in  negative  z directions  [14]. 

It  should  be  noted  that  Collin  [14]  uses  superscripts  (+)  or 
(-)  to  denote  the  direction  of  the  wave  propagation.  In  the  present 
report,  letter  a denotes  the  con^lex  amplitude  of  the  wave  traveling 
in  the  positive  z direction,  and  letter  b is  used  for  propagation 
in  negative  z direction.  These  are  the  familiar  scattering  param- 
eters [15] . 
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SECTION  II 

i’RAVELING  mWE  FORMULATION 

Lossless  nulticonductor  transmission  lines  OITL)  with  unequal 
conductors  and  inhomogeneous  dielectrics  give  rise  to  multivelocity 
quasi -TEM  waves.  As  shown  in  the  Appendix  C,  the  voltages  and  currents 
on  a MTL  are  described  by 


|V(z) 


N -j6-z  j3.z 

> = I (a.e  ^ + b.e  , (15) 

n=l  ^ ^ ^ 


N -jB.z  j3.z 
|l(z)>  = I (a.e  - b.e  )|v».> 
n=l  ^ ^ ^ 


(16) 


|<(»^>  and  are  normalized  voltage  and  current  eigenvectors,  and 

a^'s  and  b^'s  are  scattering  amplitudes  of  the  waves  traveling  in  positive 
and  negative  z direction,  respectively.  It  is  assumed  that  on  an  N 
conductor  MIL  there  are  N normal  modes,  each  of  them  having  a -distinct 
propagation  constant  3^^. 

The  notation  from  (15)  and  (16)  may  be  made  more  ccmpact  by  introducing 
the  vectors  of  incident  and  reflected  normal-mode  amplitudes 


a>  = 


^1 

f 


|b>  = 


h 


(17) 


and  by  introducing  the  diagonal  matrix  E(z)  containing  the  exponential 
functions 


jBiZ 


E(z)  » diag.  (e,  1 e ^ ) , 


(18) 
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In  the  new  notation,  (IS)  and  (16)  become; 


lV(z)>  = ^(E*(z)|a>  + E(z)|b>)  , 

(19) 

|l(z)>  = Mj(E*(z)|a>  - E(z)lb>)  , 

(20) 

where  * denotes  a complex  conjugate  number,  and  ^ and  Mj  are  matrices 
consisting  of  voltage  and  current  eigenvectors,  as  defined  by  (C-42)and  (C-43). 
(19)  and  (20)may  be  now  solved  for  |a>  and  |b>: 

E*(z)|a>  = |(MjlV(z)>  + ^|I(z)>)  , (21) 

E(z)lb>  = |(MjlV(z)>  - ^ll(z)>)  . (22) 


|a>  is  a constant  vector  containing  the  anqjlitudes  of  the  individual 

modes  as  its  components.  For  instance,a^  is  the  complex  amplitude  such  that  the 

total  power  of  the  i-th  mode  transmitted  in  positive  z direction  is 


The  entire  power  of  all  the  modes  traveling  in  the  positive  z direction  is 


• <a  a> 


(23) 


^ere<a|  denotes  a transpose  conjugate  of  |a>.  Similarly,  the  entire 
negative- traveling  power  is 

P'  = ^<b  ! b >.  (24) 

The  net  power  is  the  difference  of  the  two.  On  a uniform  MTL  there  is 
no  exchange  of  power  between  different  modes.  Each  mode  travels  with 


f 


constant  magnitude,  while  its  phase  grows  linearly  with  distance. 
The  z dependence  of  scattering  amplitudes  may  be  expressed  as 


and 


|a(z)>  * E*Cz) |a> 


|b(z)>  = E^(z)|b> 


(25) 


(26) 


where  la>  and  |b>  are  vectors  consisting  of  ccxnplex  constants.  Thus, 

the  i-th  conqxinent  of  the  vector  equation  (.25)  is 

-j  SjZ 

a^(z)  = e a^  , 

and  the  corresponding  i-th  con^wnent  of  (26)  is 

b^(z)  = e b^  . 

The  signal  flow  graph  [l6]of  the  MTL  section  of  length  Z is  shown  in  Fig.  4. 

The  i-th  mode  has  two  variables  a^(0)  and  bj^(O)  at  z=0.  Similarly,  at 

z=i  the  two  variables  of  the  i-th  mode  are  a.(£.)  and  b.(£.).  The  coefficients 

of  matrix  E* (i)  equal  to  e . Thus , a . (£)  is  obtained  by  multiplying 

-j6i£ 

a^(0)  by  the  coefficient  e . The  signal  flow  graph 

in  Fig.  4 is  extremely  simple,  since  along  the  transmission  line  there 
is  no  cross  coupling  between  different  modes. 

Figure  5 shows  a MTL  section  with  3 conductors  above  the  ground  plane. 

At  z=0,  a source  is  coupled  to  the  NTTL,  inducing  the  waves  traveling  toward 
+z  and  -z  directions.  The  amplitudes  of  the  i-th  mode  produced  by  the 
source  are  a^^  and  b^^.  The  direction  of  propagation  of  these 

waves  is  indicated  by  a wavy  arrow.  The  wave  traveling  in  +z  direction 
arrives  at  z*£^  shifted  in  phase  for  e . There,  a three -port  net- 
work terminates  the  MTL.  The  impedance  matrix  of  this  network  is 
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The  amount  of  reflection  on  the  MTL  can  be  computed  from  ^ as  follows. 
At  z»0,  the  total  amplitude  a^  of  the  i-th  mode  traveling  in  the 
positive  z direction  consists  of  two  parts:  the  source  part  a^^  and 
the  part  which  arrives  reflected  from  the  left-hand  temination  at 
z*-ly  The  total  amplitudes  of  all  waves  at  z=0  are  arranged  in  vector 
|a>.  Then,  the  vector  of  all  waves  arriving  at  z»4^  is 

|a(V>  * E*(il^)|a> 

From  (21),  this  is  furthermore  equal  to 

From  (22),  the  reflected  mode  vector  at  the  end  of  line  is 


|b()l4)  = ^i|V(44)>-  ^|l(il4)>), 


The  currents  and  voltages  at  the  load  network  terminals  are  related 


through  7^ 


|V(V>  = 


By  eliminating  |V(Jl^)>  and  ll()l^)>  from  the  last  three  equations  one 


obtains 


|b(i4)>  * |a(^4)> 


where  the  reflection  matrix  at  port  4 is 


§4  - (Mj  ^4Mj-  y)(Mj  U)- 


c 


When  the  load  impedance  matrix  is  equal  to  the  characteristic 
impedance  7^*  2^,  one  has  a reflectionless  termination,  as  can  be  verified 
by  (C-46)  and  fC-SO).  Otherwise,  is  a matrix  which  has  usually  non- 
vanishing off-diagonal  elements.  This  means  that  a single  incoming  mode 
a^  produces  a multitude  of  reflected  modes  bj(j=l,2,3,  etc.).  Thus, 
energy  transfers  from  any  one  mode  into  all  other  modes. 

This  fact  is  illustrated  in  the  signal -flow  diagram  on  Fig.  6.  On 
top  of  the  figure,  a^^^  represents  the  amount  of  energy  coupled  from  an 
external  source  to  mode  1,  traveling  in  +z  direction.  The  total  wave 
amplitude  a^  at  the  origin,  consists  of  a^^  plus  the  wave  which  was 
reflected  from  the  termination  at  z=-jl,; 


This  wave  arrives  at  port  4 as 


ajU^)  = a^e 

A portion  of  this  wave  goes  back  as  a reflected  wave  of  the  mode  1 : 

^4,11®1^V 

vdiere  is  the  coefficient  (1,1)  of  the  matrix  S^. 

Another  portion  of  the  wave  a^^Ct^)  goes  to  the  enhancement  of  the 


reflected  mode  2: 


^4, 21^1 ^^4^ 


and  so  on,  as  shown  in  Fig.  6.  For  example,  the  total  reflected  wave 
of  the  mode  2 is 

“ ^4, 21^1  * ^4, 22^2 ^4, 23^3 * 


Similar  situation  occurs  at  port  3,  where  the  load  network  Zj^is  attached. 


There,  the  reflected  wave  is 


|a(-M>  - SJb(-t,)>  , 


H*f!^^i_iik>^>‘J<**"**^*'****  •***»"  <j«>«<i'i»  If  Vi iiiiiiiiJitfiWfcjM 


!• 

t with 

'} 

I §J- .g-l  (33) 

' l; 

' The  above  discussion  was  based  entirely  on  the  frequency-domain 

considerations,  where  each  mode  was  a steady-state  sinusoidal  function 
- of  time.  Assuming  the  quasi-TBl  waves  are  non-dispersive,  an  arbitrary 

jJ  waveform  is  transmitted  by  each  mode  without  a distortion.  Therefore, 

for  mode  i,  the  wave  traveling  in  the  +z  direction  at  the  port 

4 is  just  a delayed  waveform  which  started  at  the  origin  as  a^Q(t): 


t 


ai4(t)  = a.jjCt 


(34) 


Since  the  waves  are  now  functions  of  both  time  and  position,  the  second  j 

subscript  is  used  to  specify  the  position:  0 for  origin  , 4 for  z=i^  etc,  a 

The  time-table  of  the  outgoing  and  reflected  waves  on  a 3 conductor  j 

line  is  presented  in  Fig.  7,  At  z=0,  the  three  waves  start  to  travel  to  I 

the  right:  1 

aio(t)  , a^pCt)  , a^pCt)  | 

and  similarly  the  three  waves  start  to  travel  to  the  left  | 


Each  wave  travels  with  its  mode  velocity  v^.  At  t=tj^,  the  first  wave 
arrives  at  port  4:  , 
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Its  shape  is  a shifted  shape  of  the  wave  \diich  started  from  the  origin, 
according  to  (34).  This  wave  produces  three  reflected  waves,  which  then 
travel  with  their  corresponding  velocities  back  in  the  negative  z direction. 
When  the  load  network  contains  inductances  and  capacitances,  the  shape 
of  the  reflected  waves  will  be  different  from  the  incaning  waveshape. 

For  simplicity,  consider  the  case  when  the  load  network  is  purely  re- 
sistive. In  that  case,  the  reflected  waves  will  be  of  the  same  shape 
as  the  incoming  waves.  The  amount  of  reflection  is  specified  by  (30) 
and  (31),  where  Z^^is  a purely  resistive  network,  thus  ^ has  also 
purely  real  elements.  The  reflections  from  the  left-hand  termination 
are  handled  in  an  analogous  manner. 


As  demonstrated  above  the  scattering  repre- 

sentation in  terms  of  vectors  |a(z)>  and  |b(z)>  is  very  convenient 
for  the  treatment  of  propagation  and  reflection  of  waves  on  a ^f^L. 
Nevertheless,  often  it  is  necessary  to  compute  explicitly  the  voltages 
and  currents  on  individual  conductors,  They  are  specified  by  vectors 
|V(z)>  and  |l(z)>.  The  transformation  from  scattering  variables  to 
the  currents  and  voltages  is  obtained  from  (19) , (20) , (25) , and  (26) 
as  follows 

|V(z)>  = I^(la(z)>  + |b(z)>)  , (35) 

|l(z)>  = Mj(la(z)>  - |b(z)>)  . (36) 

The  inverse  transformation  is  obtained  from  (21) , (22) , (25) , and 

(26)  : 

|a(z)>  = I (Mj|V(z)>  + ^ll(z)>)  , (37) 

|b(z)>  = i (Mj|V(z)>  - ^|l(z)>)  . (38) 

Note  that  the  components  of  scattering  vectors  |a(z)>  and  |b(z)> 
correspond  to  the  individual  normal  modes,  while  the  components  of  the 
voltage  and  current  vectors  |V(z)>  and  |l(z)>  correspond  to  the  indi- 
vidual conductors.  As  an  example,  a^  is  the  complex  amplitude  of  the 
i-th  mode.  On  the  other  hand,  is  the  voltage  of  the  i-th  conductor 
w.r.  to  ground. 
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Fig.  8. (a)  shows  a section  of  OTL  of  the  length  I,  The  voltages 


V^(z)  are  specified  between  each  conductor  and  the  ground,  where 
the  conductor  side  is  considered  as  positive.  The  currents  I^(z)  on 
the  individual  conductors  are  specified  positive  when  flowing  in 
(+2)  direction.  This  apparently  trivial  fact  is  pointed  out  because 
by  using  this  convention,  the  current  at  z=SL  points  out  of  the  MTL 
section,  which  is  not  customary  in  the  network  theory.  However,  the 
present  convention  is  found  to  be  more  natural  for  the  matrix 
manipulations  which  follow.  The  consequence  of  this  convention  is 
that  the  total  power 

P = y Re  <V|I> 

is  not  always  pointed  into  the  network,  as  customary  in  network  theory, 

but  here  P represents  the  power  flow  in  the  positive  z direction. 

Thus,  at  z=0,  P is  into  the  NTTL,  at  2=1,  P is  out  of  the  NfTL. 

The  use  of  |a>  and  (b>  variables  in  place  of  |V>  and  |I> 

variables  is  indicated  in  Fig.  8.(b).  At  each  end  of  the  NfTL  there 

is  a transforming  network  N.p,  which  transforms  the  variables  according 

to  C35)-(38).  At  z=0,  the  left-hand  terminals  of  Mp  are  the  actual 

f'fTL  conductors.  Here,  the  variables  are  |V(0)>  and  |l(0)>  . The 

other  side  of  N,p  network  has  the  mode  variables  |a(0)>  and  |b(0)>. 

Each  mode  is  represented  by  a fictitious  single  transmission  line, 

shielded  from  all  the  other  lines.  Each  single  transmission  line 

passes  the  wave  through  by  simply  adding  the  phase  shift  as 

follows  -jB-Jl 

apU)  * ap(0)e  ^ 

At  the  other  end  of  the  ^^^L,  at  z=l,  variables  |a()l)>  and  |b(ll)>  are 
transformed  back  to  variables  (V()l)>  and  ( !(!.)>  in  another  identical 
transforming  network  Np. 
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SECTION  III 

EQUIVALEOT  SOURCES  FOR  APERTURE  EXCITATION 


Each  nonnal  mode  is  described  by  its  scattering  amplitude  a^ 
(waves  traveling  in  +Z  direction)  or  b^  (waves  traveling  in  -z 
direction).  The  detailed  distribution  of  the  electric  and  magnetic 
field,  for  the  i-th  mode  traveling  in  +z  direction,  is 


* 'j®i^  - 

E^(x,y,z)  = a^e  ej(x,y)  , 


(39) 


'jB-z  ^ 

fii(x,y,z)  = a^e  ^ h^(x,y) 


(40) 


e^  and  are  the  normalized  modal  field  distributions  over  the  cross 


section  of  the  MTL.  The  total  power  transmitted  by  the  i-th  mode 
traveling  in  +z  direction  is  obtained  by  integrating  the  Poynting 
vector  over  the  cross  section: 

^i  ' <i^i(x,y)]*  ds  . 
cross  section 


For  normalized  modal  field  distributions 

[e.(x,y)  X ii*(x,y)]*ds  = 6^^  , 


If: 


(41) 


cross  section 
and  the  power  transmitted  by  the  i-th  mode  is 


pI  • 


in  accordance  with  circuit  theory. 

Fig. 9.  (a)  shows  the  junction  representing  an  aperture  on  a NITL. 

The  aperture  is  excited  by  an  external  field  from  below  the  ground  plane. 
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This  excitation  is  replaced  by  an  equivalent  pair  of  dipoles  and 

c . In  what  follows,  an  equivalent  circuit  will  be  established  for 
ey 

the  junction  between  the  two  infinitely  close  planes,  denoted  by 
L and  R (letters  stand  for  "left"  and  "right") . 

According  to  (13) » the  dipole  pair  and  c^^  excites  the  i-th 
mode  traveling  in  +z  direction  as  follows 


Si  “ I '^mx^xi^S’®^  * Sy®yi^^0'°^^  ' 


(42) 


As  indicated  previously,  a convention  in  the  present  report  is  that 
the  waves  traveling  in  +z  direction  are  denoted  by  a^^,  while  the  waves 
traveling  in  -z  direction  are  denoted  by  b^^.  Therefore,  (42)  gives  a 
source  of  the  traveling  mode  in  +z  direction,  and  is  denoted  by 
A traveling  wave  source  is  a three -terminal  device,  as  explained  in 
Appendt:^  B. 


Another  word  about  the  notation  in  (42):  x-conponent 

of  the  magnetic-field  modal  distribution  of  the  i-th  mode^  evaluated 
at  the  point  x=Xq  and  y=0,  Similarly  ey,^(xQ,0)  is  the  y-component  of 
the  normalized  electric  modal  field. 

For  the  traveling -wave  source  of  the  i-th  mode  propagating  in 
-z  direction,  the  following  is  obtained  from  (9) : 


'’si  ' 1 ' 


(43) 


Fig.  9. (b)  shows  the  modal  sources  between  the  left  and  right 
reference  planes.  This  equivalent  circuit  is  appropriate  for  an  analysis 
in  terms  of  scattering  coefficients.  The  equivalent  circuit  in  terms 


of  voltages  and  currents  can  be  now  obtained  by  a simple  matrix 
manipulation.  First,  define  the  vectors  containing  the  scattering 
sources  as  follows; 


a > = 
s 

• 

lb5>  = 

• 

t 

• 

^sN 

■ ^sN- 

(44) 


The  traveling  wave  vector,  |aj^>  traveling  in  +z  direction  out  of  the 
right-hand  reference  plane  is  a sum  of  the  source  vector  |a^>  and 
the  wave  vector  |aj^>  incoming  from  the  left  upon  the  left-hand  reference 
plane: 

I V " I V I V • 

Similarly,  the  waves  traveling  in  the  -z  direction  are  related  as 

IV  = I V IV  • (46) 

The  voltage  and  current  variables  of  the  source  junction  are 
shown  in  Fig.  9.(c).  The  voltages  and  currents  at  the  left-hand 
plane  are  defined  by  |Vj^>  and  |Ij^>  , while  the  right-hand  variables 
are  |Vj^>  and  The  current  sources  from  Fig,  9,(c).  constitute 


vector  llg>  while  the  voltage  sources  make  lVg>,  The  Kirchhoff  laws 
require 


IV* 

IV  ■ IV  ’ 

(42) 

and 

I'l** 

|I,>  = |Ir>  . 

(48) 

To  change  from  scattering  representation  to  voltage  representation, 
use  (35) 

l.v  ■ IV-IV  ' tvdv-lv  • • 
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By  virtue  of  (45)  and  (46) : 


IV  * t?v^lV‘lV^ 


The  components  of  |a  > and  |b  > are  given  by  (28)  and  (29).  The 


difference  term  is 


a . - b . = c h . (x»,0) 
SI  SI  mx  XI  ^ 0 


and  it  depends  only  on  the  magnetic  dipole  excitation.  Thus,  the 
voltage  sources  in  Fig.  9(c)  are  to  be  canputed  as  follows: 


IV  ‘ M, 


The  current  sources  are  similarly  found  as  follows 


|l5>  = Mj  (|ag>+  |bg>)  . C53, 

The  sum  term  depends  entirely  on  the  electric  dipole  moment 

* '>51  ■ 

SO  that  the  formula  for  canputation  of  the  current  sources  becomes 


1 1 > * M, 
' s —I 


“eyV'^'o-'" 


c e.,(Xrt,0) 
ey  yN^  0’  ■'  j 


In  the  above  derivation,  the  energy  storage  at  the  junction  was  not 
taken  in  the  account.  The  circuits  from  Fig.  9.(b)  and  (c)  are  thus  the 
zeroth-order  equivalents,  such  as  the  zeroth-order  equivalent  of  a single 
transmission  line  from  reference  [13]. 
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SECTION  IV 


EXAMPLES  OF  VOLTAGE  WAVEFORMS 


The  voltages  induced  on  a 2-conductor  transmission  line 
filled  with  inhomogeneous  dielectric  will  be  now  confuted  in  order 
to  illustrate  the  lise  of  the  theory  developed  thus  far.  The  system 
is  shown  in  Fig.  10.  A small  circular  aperture  of  diameter  d = 2 cm 
is  located  at  z = 0.  The  two-conductor  transmission  line  is  located 
between  z = = -7m  and  z * = 5 m. 

The  cross  section  of  the  transmission  line  is  shown  in  Fig.  11: 
it  consists  of  two  strip  conductors  of  width  w,  placed  between  the 
three  layers  of  dielectrics  denoted  and  e^.  The  dielectric 

thicknesses  are  denoted  by  h^,  h2,  and  h^.  This  parallel -plate 
model  of  the  transmission  line  is  selected  because  of  its  simplicity, 
and  it  will  be  used  to  illustrate  the  procedure  of  computing 
voltages  induced  by  an  EMP  wave.  The  electrostatic  field  within  the 
transmission  line  frcm  Fig.  11  may  be  produced  in  two  independent 
ways.  In  the  excitation  A,  a potential  V^  is  applied  to  conductor 
a,  while  conductor  b and  the  shield  are  held  at  zero  potential.  Then, 
the  fields  in  the  three  regions  are 


lA  y ^ 

In  the  excitation  B,  a potential  Vj^  is 
conductor  a and  the  shield  are  held  at 
fields  are 


applied  to  conductor  b,  while 
zero  potential.  The  corresponding 
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t t t 0 t -►  ’'b 

^IB  * ° ’ ^2B  “ FJ  ' ^3B  ■ ^y  HJ  ' 

Under  arbitrary  excitation  the  fields  in  the  individual  regions 
(i  = 1,  2,  3)  are  obtained  by  superposition: 


= V„  t. 
a la 


V ^ 
''b  ^ib 


(56) 


In  order  to  confute  the  coefficients  of  electric  induction 

matrix,  one  must  find  the  charge  per  unit  length  of  the  conductor  a. 

This  is  accoii5)lished  by  integrating  the  electric  flux  through  the 

closed  surface  S around  the  conductor  a: 
a 


where  e and  ^ depend  on  the  point  of  integration.  Using  (56) , the 
chaTg®  Qo  is  found  to  consist  of  two  parts,  one  proportional  to  V, 

ci  ci 

and  the  other  proportional  to  V^. 


(57) 


Qa  = Vaa  ^ Vab 

The  constants  of  proportionality  are  called  induction  coefficients, 
denoted  K and  K . . They  depend  only  on  the  geometry  of  the  system. 

act  aD 

For  example,  K is  confuted  as  follows: 

ad 


1 


j 


rs- — - 
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K = e, 
aa  1 


®1A  • - ^3 


2_.  • a dx 
3A  y 


x=-w/2 


K = w ^ ^ 

aa  hj 


x=-w/2 


Similarly,  the  other  induction  coefficients  are  found  as  follows: 


■'ab  = = -'•  hf 


e-  e. 


[h^  h^ 


These  coefficients  form  induction  coefficient  matrix  K 


•K  K , 

aa  ab 

"ab  ^b 


As  an  exanple  of  a symmetric  system,  the  following  dimensions  have  been 
selected: 

h,  * 2 cm,  h-  = 2 cm,  h-  = 1 cm,  w = 10  cm,  e,  = 1.0, 


e-  ' 1.0,  = 2.0,  d = 2 cm. 


The  corresponding  matrix  K is 
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■ 25  -20  ' 

= ^0  -20  25 

V . 

The  induction  coefficient  matrix  with  only  a}.r  as  dielectric  will  be 
denoted  K' ; 

' 15  -10' 

V*  s c 

- ^0  [-10  15 

By  using  matrices  K and  K’  one  can  find  the  modal  velocities,  voltage 
and  current  eigenvectors,  and  the  iji5>edance  matrix  by  the  procedure 
described  in  Appendix  C.  The  inverse  of  the  induction  matrix  L is 
directly  proportional  to  K'  as  follows 

L'^  = z - K'  (62) 


where  c is  the  velocity  of  the  light  in  vacuum.  Next,  the  eigenvectors 
and  eigenvalues  of  L’^  are  found,  and  an  auxiliary  matrix  B is  formed 
according  to  (C-15).  When  the  eigenvalues  of  B are  computed,  the  modal 
velocities  v^  are  found  as  their  inverse  square  roots,  according  to 
(C-16).  The  eigenvectors  of  B are  then  used  to  form  the  modal  matrices 
and  Mj,  according  to  (C-42)  and  (C-43) : 


•^al 

^a2 

■-6. 136  - 

•2.369' 

^1 

-6.136 

2.369 

"''al 

-0.081489 

-0.21106' 

>bl 

%2^ 

-0.081489 

0.21106 

(63) 


(64) 


Subscripts  a and  b denote  the  conductors,  and  subscripts  1 and  2 
denote  the  modes. 

Next,  modal  functions  e and  h vdll  be  evaluated.  The  aper- 

y ^ 

ture  is  placed  at  the  center  of  the  bottom  shield  conductor,  as  shown 
in  Fig.  11.  By  definition,  the  electric  field  of  the  n^^  mode  trav- 
eling in  positive  z direction  is 

^j^(x,y,z)  = a^e^(x,y)e  ^ ^ (65) 

By  selecting  a^  = 1,  the  electric  field  at  z = 0 becomes 


= e^(x,y) 


(66) 


When  = 1 3nd  b^^  = 0,  the  voltage  vector  is  obtained  from  (35) 
as  follows 


n 


■^an 

^bn 


(67) 


Thus,  in  order  to  find  the  modal  function  e of  the  mode  n (here, 

n 

n = 1 or  2) , the  potentials  on  the  two  conductors  must  be  selected 
equal  to  (J>^  and  4)^^  as  shown  in  Fig.  12.  Then  the  modal  function  is 
equal  to  the  electric  field,  according  to  (66) . Since  the  field  has 
only  the  y component,  the  result  is 


®yn  ~ *^fc.i^y2A  * ‘^bn^y2B 


(68) 


The  modal  function  ^^(x,y)  is  equal  to  the  magnetic  field  f5^(x,y,0) 
inside  the  transmission  line,  when  conductor  a carries  a current 
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and  conductor  b carries  current  The  situation  is 


shown 


in  Fig.  13.  The  currents  are  assumed  to  be  uniformly  distributed  over 


the  conductor  surfaces.  An  elementary  ccmputation  gives  the  following 
value  for  the  magnetic  field  modal  function: 


. *an''l  * 


w(h^  ^2  + h-^) 


In  the  example  treated  here,  a time  domain  response  will  be  evalu- 
ated, while  most  of  the  theory  presented  until  now  has  been  formulated 


in  the  frequency  domain.  In  order  to  use  formulas  (8)  and  (9)  for  a 


general  time  A^riation,  they  are  rewritten  as  follows: 


c = -uot 
mx  ^ m 


jOxt  _ pjint 

xs  xs 


c = EO  4 
ey  e 3t  [ ys  ys 


The  polarizabilities  and  are  given  by  (10) . 


The  incident  wave  orientation  is  specified  by  angles  a and  0 as 


shown  in  Fig.  3.  Then  for  a TM  polarization,  the  field  coiiponents  of 


inportance  are 


= 2E„F(t)sin9  , 
ys  0 ^ ’ xs 


F(t)sina 


An  arbitrary  time  variation  of  the  incident  plane  wave  is  described  by 


function  F(t).  For  EMP  wave,  a simple  function  is  selected  as  follows 


-a.t  -B.t 
F(t)  . e ^ - e ^ 


****.•*•  1 * 

* » . 1 ■ . 

'ir/.;  *.  * •• 

71 ..  ‘..'w 

^':+  + 

^a\ 

+:  d-  ■+•■+  ■+  + 

•e- 

or 

0 

Fig.  12.  Evaluation  of  the  electric  field  on  the  line  from  Fig.  11 
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Fig.  13.  Evaluation  of  the  magnetic  field  on  the  line  from  Fig.  11 


I 

having  the  following  derivative 

-a,t  -B-,t 

F'(t)  = -a^e  ^ + 6^e  ^ . (74) 

The  following  parameters  of  the  incident  plane  wave  have  been 
selected  for  the  present  example; 

Eq  = lookV/m,  a = 30“,  0 = 45°,  = 3.10^s‘^,  3^  = 10®s‘^. 

The  electric  and  magnetic  current  moments  are  now  computed  by  (70) 
and  (71) : 

'^5) 

Cgy=  2 e^EQF'(t)  sine  (76) 

These  values  are  to  be  substituted  in  (42)  and  (43)  in  order  to  compute 
the  traveling  wave  sources  for  different  modes  on- the  multiconductor 
line.  For  the  present  exanple,  the  modal  functions  e^.^^  and  h^^  which 
are  needed  in  (42)  and  (43)  are  determined  by  (68)  and  (69) . The  ini- 
tial anplitudes  of  the  traveling  ^vave  sources  can  be  now  computed  and 
the  multipole  reflection  traced  down  in  a manner  indicated  in  Fig.  7. 

As  long  as  the  loading  at  each  end  of  the  multiconductor  line  consists 
of  pure  resistances,  such  as  in  Fig.  10,  the  wave  shape  remains  unchanged 
after  each  reflection,  and  is  specified  by  the  function  F'(t)  from  (74). 
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In  the  first  example  to  be  computed,  the  following  loading 


resistances  have  been  selected: 

\l  “ '^aR  “ ‘HiL  ' \r  “ ^cL  “ \r  “ 

The  load  impedance  matrix  of  a T network  is  then  obtained  as 


\ll  = ^aL  ^ 


^L12  “ ^cL  ■ 


^L22  ' ^cL  ' 


The  corresponding  scattering  matrix  is  obtained  from  (33) 


' 0.43488  0.17495' 

-3  " 

0.17495  0.90584^ 

Since  the  load  resistances  on  the  right-hand  side  are  the  same  as  on 
the  left-hand  side,  = S^.  For  each  incident  wave  a^  at  the  port  4, 
the  reflected  waves  are  computed  by 


b.  = S. .a.  for  j = 1,2 

j :i  X ’ 


Then,  the  voltages  are  obtained  from  (35) 


1V>  = ^ (la>  + lb>) 
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The  velocities  of  the  two  waves  are  = 3*10  m/s  and 

g 

^2  “ 2.236  • 10  m/s.  Fig.  14  shows  the  voltages  on  the  right- 
hand  end  of  the  transmission  line.  is  the  voltage  on  the  con- 
ductor a , while  V2  is  the  voltage  on  the  conductor  b.  First 
pulse  arrives  at  t^  = = 16.6  ns.  Note  that  this  is  even 

mode,  since  the  polarity  of  the  pulse  is  the  same  on  both  conductors. 
Shortly  afterwards,,  at  t2  = ^^4/^2  = 22.4  ns,  the  odd  mode  arrives 
producing  a positive  pulse  of  and  a negative  pulse  of  V2.  The 
next  arrival  is  the  group  of  four  waves  which  are  reflected  from 
the  left-hand  end  of  the  line,  and  the  process  is  continued  through 
multiple  reflections  bouncing  back  and  forth  on  the  line. 

At  each  of  the  bounces,  some  of  the  energy  is  lost  in  resis- 
tances terminating  the  line,  so  the  process  gradually  dies  off  as 
seen  in  Figure  14.  If  one  end  of  the  line  is  terminated  in  a resis- 
tance matrix  equal  to  the  characteristic  inpedance  matrix  there 
are  no  reflections  from  that  end.  In  the  next  exanple,  we  terminate 
the  right-hand  end  of  the  NfTL  by  a matched  inpedance.  The  necessary 
resistances  are  obtained  from  ^ as  follows.  First  ^ is  computed 
by  (C50) 

43.260  32.036  ' 

32.036  43.260  , 

Then,  the  resistances  of  the  T-network  are  obtained  as 

*^aR  * ^011  ' ^012  “ 11*2240 
\r  " ^022  ■ ^012  “ 11*2240 
^cR  “ ^012  ^ 32.0j6O 
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8b.OO  1^0.00  240.00  320.00  400.00 

TIME  IN  NRNOSECONDS 


sb.OO  1^0.00  240.00  iEoToO  7oO .00 

TIME  IN  NRNOSECONDS 


Voltage  waveforms  on  a 2 -conductor  line:  moderate  mismatch. 

i2 


A 


t 


With  such  a matched  termination  on  the  right-hand  side,  the  resulting 
voltages  take  the  shape  as  shown  in  Fig.  15.  The  first  arrival  of 
the  odd  and  even  mode  is  similar  to  the  previous  situation  in  Figure  14. 
Also,  the  next  four  pulses  are  similar  to  those  from  Figure  14,  but 
after  that  there  are  no  more  waves  left  on  the  line. 

Another  extreme  situation  occurs  when  there  is  no  attenuation  , 
because  the  terminations  consist  of  open-circuits  or  short  circuits. 

In  Figure  16,  the  following  terminations  have  been  selected 

\l  ^ *^aR  \l  " \r  " ^cL  “ ^cR  ^ ^ 

There  is  little  attenuation  of  the  voltage  waveshapes  in  the  first 
300  ns,  and  the  multiple  bouncing  on  the  transmission  line  continues 
for  a long  time. 

The  conqjuter  program  evaluates  the  individiaal  pulse  arrivals  in 
the  close  analogy  with  the  time-table  from  Figure  7.  If  the  ^f^L  can 
support  N different  modes,  each  of  those  produces  N other  reflected 
modes  at  each  bounce.  Very  soon  the  storage  requirements  become 
prohibitive.  For  a 2-conductor  line  of  the  total  length  12  m as  in 
Figure  10,  there  are  510  pulse  arrivals  within  300  ns.  If  the  number 
of  conductors  is  increased  to  three,  there  are  9840  arrivals  which  are 
to  be  stored  before  sorting  them  in  chronological  order.  It  is  obvious 
that  this  procedure  is  practical  only  for  computing  the  early  time 
responses . 

An  example  of  the  three  conductor  parallel -plate  ^^^L  is  shown 
in  Figure  17.  The  voltages  on  the  individual  conductors  are  denoted 
by  V^,  V2,  and  V3. 
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Fig.  17.  Voltage 


All  the  exanples  computed  in  this  Section  utilize  a parallel- 
plate  model  of  the  MTL  with  inhomogeneous  dielectric,  because  of 
the  availability  of  simple  formulas  for  the  evaluation  of  induction 
coefficient  matrix  K and  of  the  normalized  modal  functions  ey,(xQ,0) 
and  h^Cxp.O).  However,  the  same  procedure  of  computing  the  voltage 
waveforms  may  be  applied  to  any  other  set  of  data  for  K,  e^,  and 
vdiich  may  be  obtained  by  a numerical  solution  of  the  arbitrary 
shapes  of  conductors.  Several  numerical  procedures  for  ccxnputation 
of  the  induction  coefficient  matrix  K have  appeared  in  the  recent 
literature  (see  [16]  to  [18]).  If  these  methods  are  supplemented  by 
computation  of  OyCx^.O)  and  h^(xQ,0)  at  the  aperture  center,  the  rest 
of  the  computations  of  voltage  waveforms  described  in  this  Section 
is  applicable  to  arbitrary  conductor  shapes  and  sizes. 


16.  W.  T.  Weeks,  "Calculation  of  Coefficients  of  Capacitance  of 
Milticonductor  Transmission  Lines  in  the  Presence  of  Dielectric 
Interface,"  IEEE  Transactions  Microwave  Theory  Techn.  Vol.  MTT-IS, 
pp.  35-43,  Jan.  1970. 

17.  J.  C.  Clements,  C.  R.  Paul,  A.  T.  Adams,  "Confutation  of  the 
Capacitance  Matrix  for  Systans  of  Dielectric-Coated  Cylindrical 
Conductors,"  IEEE  Transactions  on  Electrcanag.  Ccwnpat.  Vol.  EMC-17, 
No.  4,  pp.  238-248,  Nov.  1975.  Also  see  Correction  in  EMC-18 

No.  2,  pp.  88-89,  May  1976. 

18.  C.  R.  Paul,  "Computation  of  the  Transmission  Line  Inductance 
and  Capacitance  Matrices  from  the  Generalized  Capacitance 
Matrix,"  IEEE  Transactions  Electromag.  Conpat.,  Vol.  BC-18. 

No.  4,  pp.  175-183,  Nov.  1976. 
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SECTION  V 


FIRST-ORDER  EQUIVALENT  CIRCUIT  OF  THE  SMALL  APERTURE 

In  Fig.  13  , a set  of  incident  waves  coining  from  the  left  is  described 
by  aperture  region  is  located  between  two  planes  denoted  L 

and  R.  There  is  no  incident  wave  coming  from  the  right,  |b„>  = 0.  The 
waves  |aj^>  excite  the  pair  of  dipoles  c^^,  and  c^.  The  excitation  field 
produced  by  the  j-th  incident  mode  is 


(77) 

^3  “ 

(78) 

This  is  an  internal  field,  according  to  Fig.  2.  From  (8)  and 

the  dipole  moments  produced  by  the  j-th  mode  are 

(9), 

(79) 

Syj  ' 

(80) 

Summed  over  all  the 

incident  modes: 

N 

(81) 

(82) 

where  argument  (Xq,0)  has  been  ommitted  for  brevity. 


These  dipoles  produce  the  outgoing  waves  • Their  i-th  component 
is,  according  to  (13)  and  (14): 


T 


W'liJViWOiii 


(92) 


R = j H + j E , 

T = U - j H + j E . (93) 

The  signal  flow  graph  [16]  is  in  Fig.  19.  Also  shown  are 
sources  |ag>  and  |bg>  from  the  zeroth-order  equivalent  circuit.  Figure 
19  is  a complete  first-order  equivalent  circuit  for  scattering  represen- 
tation. The  corresponding  immittance  representation  will  be  derived  next 


Fig.  20  shows  a small  parallel  perturbation  on  a NfTL.  The  perturba- 
tion is  described  by 


and 

M 

A 

(94) 

|l2>  = Y |Vj^>  . 

(95) 

Kirchhoff  current 

law  requires 

IV  = IV"  IV  • 

(96) 

Change  variables  |V>  and  |l>  into  |a>  and  |b>  according  to  (35)  and 
(36).  Then,  (94)  and  (96)  become; 

V'v  * ' V'v  * • 

Mi(|aj^>  -|bL>)  = YM^daL^  ’ 

Multiply  the  first  equation  from  the  left  by  and  the  second  equation 
by  Mj^.  Use  (C-46)  and  (C-47)  to  obtain 
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Fig.  20.  Small  parallel  admittance  perturbation  on  MTL 


liu>  II«>=IIl> 


lv.>  l'^>'  |v,> 


/m 


Fig,  21.  Small  series  impedance  perturbation  on  NITL. 
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I V IV  = iv  IV 


\a^>  - \b^>  = Y.\a^>  ^ xl  V " I V ' IV  ' 


where  the  normalized  admittance  matrix  is  introduced  as 


y = M*'  Y M 
^ - -V 


By  subtracting  the  above  two  equations  one  obtains 


\h^>  = -(2U+x)'^y|aL>  + 2(2U+y)'l|bj^> 


When  the  normalized  admittance  matrix  is  "small"  [19]  it  is  possible  to 
use  the  first  two  terms  from  Neumann's  series: 

(2U  + - IZ)  • 

Retaining  only  the  linear  terms  in  y,  (98)  is  approximately  given  by 

|bL>  ~ • yz  Iv  ■"  i>^^lV  • 

Next,  consider  a small  series  perturbation  in  Fig.  21,  described  by 


IV  = IV  ' 

|Vj>  = Z|Ij^>  , 

IV  ' lVi>  ^ IV 


(100) 

(101) 


(102) 


Expressed  by  scattering  variables 


Mj(|aj^>  - |bj^>)  = Mj(|aj^>  - |bj^>)  , 

^^Iv  IV^  “ i^idv  ■ iv^  -v^iv  iv^ 


Introduce  the  normalized  impedance  matrix 


1 “ Mi  iMi  • 


Using  the  same  approximations  as  in  (99)  one  obtains 


(103) 


IV  " iv  ^ - i^^iv  • 


(104) 


Now  add  both  series  and  parallel  perturbations  as  in  Fig.  22.  The 


small  reflections  are  simply  added  as  follows 


|1Y>  . |ai,>  * (U  - f i - 1 1)  |bR>  . 


Invoking  reciprocity  and  symmetry  of  the  junction 


h-h 


u ■ • iz 
-Zi 


Comparing  (106)  with  (91) , one  concludes 


jH  + jE  = j(z  - x)  , 
-jH  + jE  = - 1(2  + x)  • 


Solving  for  ^ and 


z = j2H 


L = -j2E 


(105) 


(106) 


(107) 

(108) 


Fig.  22.  MTL  with  both  series  and  parallel  perturbations. 


|Vs>  I 

I 


Fig.  23.  Equivalent  circuit  of  the  aperture  junction  with  sources. 
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From  (89)  and  (90)  it  can  be  seen  that  all  the  elements  of  H and  of 
E are  real  and  proportional  to  frequency  w.  Thus,  the  elements  of 
^ are  represented  by  self  and  mutual  normalized  inductances 

z = (109) 


where : 


These  inductances  may  be  positive  or  negative,  depending  on  the  signs 
of  h . and  h . . 

XI  xj 

Analogously,  the  normalized  admittance  matrix  consists  of  self- 
and  coupling  capa  Ltances: 

X = -j  0)  c (111) 


where  the  elements  of  c are 


(112) 


The  negative  sign  signifies  that  the  capacitance  is  negative  if  the  pro- 
duct ®yi®yj  comes  out  to  be  positive.  A note  of  explanation  is  necessary 
on  the  meaning  of  the  totaling  capacitance  c^j  (when  i ^ j).  Such  a cir- 
cuit element  does  not  exist  in  lumped-circuit  theory.  It  signifies  that 
the  current  at  the  port  i is  proportional  to  the  rate  of  change  of  vol- 
tage at  the  port  j,  the  constant  of  proportionality  being  defined  as  a 
coupling  capacitance. 

Finally,  to  obtain  the  actual  impedance  matrix  2,  (103)  has  to  be 
denormal ized  as  follows 


(113) 
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I 


1 


Similarly,  the  denormalized  admittance  matrix  is 
Y = Mj  X = - jo)  Mj  c ^ • 


(114) 


Note  that  all  the  matrices  on  the  right-hand  sides  of  (113)  and  (114)  are 
real,  thus  the  equivalent  circuit  is  easy  to  interpret. 

When  also  the  voltage-  and  current  sources  from  (52)  and  (55)  are 
included  in  the  circuit,  the  complete  first-order  equivalent  circuit  takes 
the  form  shoun  in  Fig.  23. 

This  representation  is  valid  below  the  first  resonant  frequency 
of  the  aperture.  For  a circular  aperture  of  radius  a the  lowest  resonance 
appears  at  [20] 

k a = 1.841 

Thus,  the  representation  is  valid  for  frequencies 

. „ 87.9*10^ 


For  an  aperture  of  radius  a = 10  cm,  the  equivalent  circuit  is  valid  for 
f«  880  Mlz.  The  validity  of  the  equivalent  circuit  could  be  further  ex- 
tended in  the  region  close  to  the  aperture  resonance  by  the  methods 
described  in  [21],  which  will  not  be  pursued  here. 

The  circuit  diagram  of  matrix  Z for  a 3-conductor  system  is  shown  in 
Fig.  24a.  The  self  inductances  are  denoted  by  and  L22,  and 


[20]  R.  E.  Collin,  Foundations  for  Microwave  Engineering,  New  York: 
McGraw-Hill,  1966,  p.  111. 

[21]  G.  L.  Matthaei,  L.  Young,  E.  M.  T.  Jones,  Microwave  Filters, 
Impedance -Matching  Networks,  and  Coupling  Structures,  New  York: 
NicGraw-Hill,  19647  P-  242. 
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the  mutual  inductances  are  denoted  by  L^2*  ^13*  ^23*  voltage 

is  then  given  by 


ju,(L^lVi  . ^ LJ3V3) 


12'2 


13’ 


Analogous  expressions  are  valid  for  V2  and  V3.  This  circuit  is  familiar 
from  the  conventional,  circuit  theory  and  it  does  not  require  further 
coiranents.  However,  the  circuit  representation  of  the  matrix  Y is  sane- 
what  unconventional.  The  self  capacitances  are  denoted  by  C22  » 
and  C23,  while  the  coupling  capacitances  are  denoted  by  C-^2>  ^23’ 

Cj^3.  The  current  is  given  by 

h - * W 

and  similar  expressions  may  ’.s;  written  for  I2  and  I3.  The  negative 
sign  signifies  that  all  these  capacitances  are  negative,  as  compared 
with  the  conventional  capacitances.  From  the  point  of  view  of  energy 
balance,  the  negative  capacitance  is  an  active  element,  capable  of 
supplying  the  energy  to  the  rest  of  the  network.  This  may  be  considered 
as  an  inconsistency  of  the  equivalent  circuit  of  the  aperture,  but  it 
will  be  shown  that  these  negative  capacitances  are  extremely  small. 
Furthermore,  the  negative  capacitance  can  be  thought  of  as  the  element 
which  represents  missing  capacitances  on  a uniform  transmission  line  which 
is  perturbed  by  the  aperture. 

The  values  of  the  elements  in  the  equivalent  circuit  from  Fig.  24 
will  be  next  confuted  for  a 2-conductor  parallel-plate  NfTL  described 
in  Section  IV.  For  a circular  aperture  of  diameter  d = 2 cm  the  electric 
modal  functions  are  confuted  from  (68) 
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eyl^O-O) 


-6.1358 


613.58 


evoCx^.O)  = = -236.90 

^ 1.10’^ 


The  values  of  the  magnetic  modal  functions  at  the  position  of  the 
aperture  are  found  from  (69) 


h^j(xQ,0)  = -0.81491  ^c2^^0’°^  " 0.42212 


From  (110)  and  (10)  the  normalized  inductances  are  obtained  as 
I = 1.1127  • 10"^^,  1^2  " 2.9855  • lO’^^  , = -5.7636  • lO'^^ 


For  the  frequency  1 GHz,  the  normalized  reactances  are 

0)^^^  = 6.9913  • 10"^  , 03)122  " 1*^758  • lO'^  , 03)lj^2  = -3.6214  • 10‘^ 

Similarly  the  normalized  susceptances  in  parallel  with  the  MTL  are 
found  from  (112)  and  (10) 


03C^^  = 1.3963  • 10‘^  , o3C^2  " -5.3909  • lO'^  , uic^z  " 2.0814  • lo'^ 


Therefore,  even  at  the  highest  frequency  of  interest  for  EMP  calculations, 
the  normalized  reactances  and  susceptances  are  small  numbers.  Thus  it 
is  expected  that  they  do  not  cause  appreciable  reflections  when  a wave 
is  propagating  along  the  NfTL. 


» 


Consider  that  £in  incident  wave  |aj^>  is  COTiing  from  the  left 
toward  the  equivalent  circuit  of  the  aperture  in  Fig.  22.  Assuming 
that  the  right-hand  side  of  the  equivalent  circuit  is  terminated  in 
an  infinitely  long  NfTL,  it  is  of  interest  to  find  |bj^>  in  terms  of 
Since  for  an  infinitely  long  MIL  the  vector  |bj^>  vanishes, 
the  equation  (106)  gives 


\h^>  = h(z-y)  |aj^> 


Therefore,  the  reflections  due  to  the  presence  of  the  equivalent  cir- 
cuit on  a uniform  ^^^L  are  determined  by  reflection  matrix  £ 


For  the  case  under  consideration,  the  reflection  matrix  is 


1=  j 


1.0478  • 10 
-9.0123  • 10 


-9.0123  • 10 
3.9572  • 10 


It  can  be  seen  that  the  reflected  amplitude  of  the  first  mode  is  only 
1.05%  of  the  incident  amplitude  for  the  same  mode.  The  other  reflections 
are  even  smaller  than  1%.  At  lower  frequencies  all  coefficients  of  the 
reflection  matrix  are  proportionally  reduced.  This  fact  justifies 
the  procedure  from  Section  IV  of  computing  the  voltage  waveforms  by 
entirely  neglecting  the  presence  of  inductances  and  capacitances  in 
the  equivalent  circuit. 
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It  is  of  interest  to  investigate  the  equivalent  circuit  of  the 
aperture  on  a single-conductor  transmission  line.  The  geometry  of  the 
problem  is  such  as  specified  by  Fig.  3:  a round  wire  of  radius  r is 
placed  in  parallel  with  the  ground  plane,  so  that  the  center  of  the 
wire  is  elevated  above  the  ground  plane  for  distance  d.  The  circular 
aperture  of  radius  a is  located  atx  = XQ,  y = 0,  z = 0. 

When  N = 1,  it  follows  from  formulas  in  the  Appendix  C that 
matrices  Mj  and  ^ reduce  to  single  scalars: 


where  C'  and  L'  are  the  distributed  capacitance  and  inductance  of  the 
single-\d.re  transmission  line.  When  the  subscripts  i and  j are  omitted, 
the  normalized  inductance  from  (110)  becomes 

and  the  normalized  negative  capacitance  from  (112)  becomes 

c = ea  e^(x-,0) 
e y^  0’ 

The  normalized  modal  functions  h and  e for  the  TENl  mode  of  a single 

X y 

wire  above  the  ground  plane  have  been  derived  in  [13]  as  follows 

[13]  D.  Kajfez,  "Excitation  of  a Terminated  TDI  Transmission  Line  Through 
a Small  Aperture,"  Interaction  Note  215,  July  1974. 


vrtiere  h denotes  the  reduced  height  of  the  wire: 


In  order  to  obtain  the  actual  values  of  the  equivalent  inductance  and 
capacitance,  the  values  are  to  be  denormalized  according  to  (113)  and 
(114): 

2 

Z = jcoL  = jcoM^  z 
Y = -jcjC  = 


The  values  of  L and  C are,  therefore 


Zq  = 60  cos  h*^  (^) 

The  similar  formulas  have  been  derived  by  Lee  and  Yeing  [22].  Their 
formulas  (20, a)  reduce  to  (115)  and  (116)  in  the  case  of  a thin  wire 
(r  « d,  h * d) . For  thick  wires  it  is  believed  tliat  (US)  and  (116) 
are  more  accurate,  because  the  exact  TEM  modal  functions  have  been 
used  in  the  derivation. 


% 


[22]  K.  S.  H.  Lee,  F.  C.  Yang,  "A  Wire  Passing  by  a Circular  Aperture 
in  an  Infinite  Ground  Plane,"  Dikewood  Corporation,  February  1977. 
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SECTION  VI 


INTERACTION  BEnvEEN  AN  APERTURE  AND  A SINGLE  WIRE 


In  this  section  we  examine  the  validity  of  the  assumption  that  the 
dipole  moments  in  the  aperture  can  be  determined  from  the  plane  wave  exciting 
the  aperture  ignoring  the  presence  of  the  NTL,  This  assumption  has  been 
used  to  derive  the  model  for  aperture  coupling  to  the  lines  in  the  preceding 
sections.  For  sinplicity,  we  treat  only  a single  wire  line  over  a ground 
plane.  It  is  possible,  in  principle,  to  extend  the  analysis  to  treat  a MTL 
backscattering  into  the  aperture  by  super inqjosing  the  backscatter  from  the 
individual  lines,  accounting  for  the  mutual  interaction  of  the  various  con- 
ductors. The  general  situation  is  too  difficult  to  treat  here,  however.  As 
an  additional  simplification,  we  assume  the  conducting  line  is  bare. 

Referring  to  Figs.  2 and  3,  in  addition  to  the  plane  wave  fields 
and  exciting  the  aperture,  we  must  now  consider  aperture  fields 

and  which  are  the  fields  reradiated  from  the  transmission  line. 
These  fields  are  linearly  related  to  the  dipole  moments  in  the  aperture  and 
we  show  in  the  following  that  they  can  be  \>n"itten  as 


H = t c 
xs  m mx 


(117) 


gint 

ys 


t c 
e ey 


(118) 


where  'tj^(tg)  is  the  x(y)  component  of  the  short  circuit  magnetic  (electric) 
field  at  the  aperture  reradiated  from  the  line  due  to  the  appropriately 
directed  unit  magnetic  (electric)  current  moment  in  the  aperture.  In  the 
absence  of  any  other  sources  in  the  interior  region,  Eqs.  (8)  and  (9) 
become  with  (117)  and  (118), 


mx 


„ext 


(119) 
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j 


I 


l+j(jjea  t 


It  is  easily  seen  that  the  coupling  from  the  line  to  the  apert 
be  neglected  if  it  can  be  shown  that 


In  the  following  we  determine  t^  and  t^^^,  the  aperture  fields  s 
the  wires  due  to  unit  electric  and  magnetic  dipoles  in  the  ape 
Beginning  with  the  dipole  sources  c^  and  c^^  (see,  e.  g., 
and  a single  wire  in  the  internal  region,  the  ground  plane  is 
that  the  dipole  and  wires  are  imaged  as  in  Fig.  25.  Note  that 
doubles  the  current  moments  of  the  dipole  sources.  Transverse 
nates  p = (p,<{))  and  p^  = (p^,4)^)  measured  from  the  wire  axis  a 
wire  axis,  respectively,  are  also  introduced  in  Fig.  25.  The 
establishes  the  location  of  the  aperture  with  respect  to  the  w 
The  fields  due  to  the  dipoles  in  the  absence  of  the  condu 
determined  from  free  space  magnetic  and  electric  vector  potent 


respectively,  where 


I 


The  currents  on  the  wire  and  its  image  are  assuned  to  be  entirely  z-di 

but  in  opposite  directions.  Hence  they  determine  a magnetic  vector  po 

A which  can  in  turn  be  used  to  determine  all  field  quantities.  Bound 
z 

conditions  require  the  z-component  of  electric  field  to  vanish  at  the 
of  the  upper  wire  in  Fig.  25;  however,  it  must  vanish  everywhere  inter 
the  wire  as  well  and  if  we  restrict  our  observation  to  the  wire  axis, 
only  the  total  or  average  current  I(z)  from  the  upper  wire  contributes 
the  vector  potential  there.  This  choice  of  the  observation  point  thus 
out  of  the  problem  any  circumferential  variation  in  the  current  on  the 
wire;  the  circumferentially  varying  currents  on  the  image  wire  do,  how 
contribute  to  the  vector  potential  at  p = d.  But  if  the  wire  is  thin  a 
sufficiently  far  from  the  ground  plane  in  terms  of  the  wire  radius,  i. 

kr  « 1, 

r « d. 


then  the  image  wire  may  essentially  be  replaced  by  a line  source  along 
image  wire  axis  which  carries  i current  I(z)  directed  opposite  to  that 
upper  wire.  In  other  words,  under  conditions  (125)  and  (126),  if  the 
wire  is  replaced  by  multipole  line  sources  located  at  p^=^,  the  c ntr 
of  the  dipole,  quadrupole  and  higher  order  multipole  terms  can  be  neg 
compared  to  that  of  the  monopole  current  filament  I (z) . Thus  the  magr 
vector  potential  along  the  upper  wire  axis  due  to  the  wire  and  image  c 
is  approximately 


Along  the  axis  of  the  upper  wire  must  vanish,  that  is 
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(128) 


a ‘f  (k^+VV-)(a  A +a  A,)  • ^x(a  F ) 1 = 0 

z L JtiJye  ^ ^ y ey  z z*  z x 

p = 0 

or  writing  out  the  desired  conqjonents , 


1 

2 1 
-it  + k2 

A 

_ 1 ey 

jcoye 

dz 

jojye  dzdy 

(129) 


Eq.  (129)  with  (123),  (124),  and  (127),  is  an  integro- differential  equation  for 

the  current  I(z)  induced  on  the  transmission  line  by  current  moments  c and 

ey 

c in  the  aperture.  In  order  to  determine  the  current  we  introduce  the 
mx  ^ 

Fourier  transform  pair 


'l^(z)  = -= 
/2tF 


t|;(k^)e  ^ dk^ 


(130) 


1 C 

ijj(z)  e dz 

^2?  , , 


(131) 


where  the  wavy  line  indicates  Fourier  transformed  quantities.  Writing  the 
current  in  (127)  as  an  inverse  transform  and  making  use  of  the  identity 


-jk/p  +(z-z')‘ 
/p^+(z-z')^ 


00 

^ J ” 


f2I  "jk  (z-z') 

J')(kpP)e  dk^ 


(132) 
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where  k = ,Rek  >0,  lmkp<0,  (127)  becomes 

r\  T A ^ 


4j(2Tr)' 


I(k,)[Hj2)(k'r) 


- H^2^C2kM)]  e 


-jCk^-k^z' 


(133) 


-jk'z 

X e dk  dk ' dz ' 
z z 


Noting  that 


03 


(k^-kpz* 


2TT5(k^-kp 


(134) 


(133)  simplifies  to  the  transform  representation 


f2'l 

- Hj^'’(2kpd)]e  ^ dk^ 


(135) 


Siiilarly,  we  express  the  right  hand  side  of  (129)  in  transform  representa- 


tion: 


3^  ^ ^ey  ® 


3z3y  3z3yj_  2Tr|?-^J  JH 


w 

tj?  any  J ^ p=tl 


pc  d 
ey 


4it/  2.  ,2 

v* 


J '''''' 


C136) 


Substituting  (135),  (136)  and  (137)  into  (12'i)  and  taking  the  Fourier 
inverse  of  the  resulting  equation  yields  finally  the  transform  of  the  current, 


(138) 


where  I (k  ) and  I (k  ) are  identified  as  the  partial  currents  arising  from 
m z 6 2 1 T • 

the  corresponding  source  terms  c^^  and  c^^,  respectively.  Note  that  is 

an  even  function  whereas  is  an  odd  function  of  the  variable  k^. 
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The  total  line  current  can  be  obtained  by  taking  the  Fourier  inverse 
of  (138).  Our  interest,  however,  is  in  obtaining  the  fields  that  the  wire 
• and  image  currents  scatter  into  the  aperture  and  these  can  be  found  most 

easily  from  the  vector  potential  expressed  in  terms  of  the  current  transform. 
The  vector  potential  at  an  arbitrary  point  is  obtained  by  again  treating 
both  the  line  and  its  image  as  current  filaments  at  p=^  and  p^=^,  respectively: 


00 


-jk(p+a  (z-2')  I 
e ^ 
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(141) 


where  we  have  used  the  even  and  odd  properties  of  and  to  appropriately 
sinplify  the  integrals  (140)  and  (141). 

Returning  to  the  examination  of  conditions  (121)  and  (122),  we  first 
note  that  the  four  parameters  k,  Xq,  d,  and  r in  the  integrals  can  be 
reduced  to  three  parameters  if  all  distances  are  measured  in  terms  of  wave- 
lengths. Thus,  we  examine 


m m 


(SF)^(k&)^kd)^t^(kxQ,  kd,  kr) 
(kxj^  ♦ (kd)^ 


(142) 


[11]  S.B.  Cohn,  "Determination  of  Aperture  Parameters  by  Electrolytic- 
Tank  Measurements",  Proc.  I.R.E.,  Vol.  39,  pp.  1416-1421,  Novanber 
1951. 

[12]  S.B.  Cohn,  "The  Electric  Polarizability  of  Apertures  of  Arbitrary 
Shape,"  Proc.  I.R.E.,  Vol.  40,  pp.  1069-1071,  September  1952. 
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(147) 


Since  our  assumption  that  the  aperture  may  be  replaced  by  current  moments 

c and  c requires  that  the  aperture  be  small,  say  less  than  a tenth  of  a 
cm  'z 

wavelength,  then  (kJl)  < .25  and  since 


(SF)o  n,  < 1 
^ e,m 


kd" 


(kxg)^  + (kd) 


I ^ 1 


(148) 

(149) 


we  need  only  demonstrate  that 

1^1  « 1.  tisw 

|i^l  « 1,  (151) 

to  show  that  the  wire  does  not  sufficiently  excite  the  aperture  that  its 
effect  needs  to  be  accounted  for  in  computing  the  aperture  dipole  moments. 

The  integrals  (146)  and  (147)  are  numerically  evaluated  by  (1)  noting 
the  syrmetry  of  the  integrals  about  a = 0;  (2)  integrating  numerically  using 
Gauss -Legendre  integration  on  the  interval  0 <.  a <.  2,  with  a singularity 
of  the  form  l/(ct-l)  removed;  (3)  adding  a contribution  from  the  branch  point 
at  a = 1;  and  (4)  integrating  numerically  on  the  interval  2 i a i “ using 
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1 
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1 


i 


Gauss- Laguerre  integration  which  assumes  an  exponentially  decaying  integrand. 

The  appropriate  rate  of  exponential  decay  in  the  integrand  can  be  determined 
from  the  large  argument  approximations  to  the  Hankel  functions  appearing  in 
(146)  and(147). 

From  the  numerical  calculations,  Figs,  26-33  show  that  over  a wide  range 
of  parameters  kx^,  kd,  and  kr,  Eqs.  (150)  and  (151)  are  indeed  satisfied.  How- 
ever, wiien  the  wire  is  directly  above  and  close  to  the  aperture  (kx^  = 0 and 
kd  small) , (150)  and  (151)  are  not  satisfied.  If  the  wire  is  close  to  the 
aperture,  however,  the  aperture  must  be  small  in  order  for  the  dipole  moment 
representation  to  remain  valid.  Calculations  by  Lee  and  Yang  [22]  indicate 
that  the  dipole  moment  representation  of  a circular  aperture  (without  considering 
backscattering  from  the  wire)  is  accurate  to  about  10?i  only  if  2kJl  < kd  and  it 
seems  reasonable  that  this  upper  bound  on  the  maximum  aperture  dimension  would 
hold  for  other  aperture  shapes  as  well.  Using  this  as  an  upper  bound  on  kJl 
and  noting  that  t^^  and  t^  are  largest  when  kx^  0,  it  is  sufficient  to  show 
that 


'kdl 

iTj 


3. 

|tm(0,kd,kr) 


« 1 


(152) 


tj^(0,kd,kr) 


« 1 


(153) 


in  order  to  exclude  consideration  of  wire-to-aperture  coupling.  These  quanti- 
ties are  plotted  in  Figs.  34-35.  The  figures  show  that  the  effect  of  the  wire 
on  the  aperture  dipole  moment  is  negligible  whenever  the  representation  of  the 
aperture  by  dipole  moments  is  valid. 

This  completes  the  task  set  for  this  section,  but  it  is  interesting  and 
instructive  to  return  to  the  current  transform  (138)  and  use  it  to  find  the 
current  on  the  line  by  inverse  transforming: 


[22]  K.  S.  H.  Lee,  F.  C.  Yang,  "A  Wire  Passing  by  a Circular  Aperture  in  an 
Infinite  Ground  Plane,"  Dikewood  Corporation,  February,  1977. 
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Fig.  26.  Magjutude  of  for  a wire  of  radius  kr  = 0.1. 


KR=  0.1000 


a — O KD= 
O — © KD= 
A — A KD= 

i h KD= 

X — X KD= 

❖ ^ KD= 

^ ^ KD= 
X — X KD= 
Z — Z KD= 


0.2000 

0.3000 

0.4000 

0.5000 

0.6000 

0.7000 

0.8000 

0.9000 

1.0000 


KR=  0.0010 


O — O KD= 
O — © KD= 
A — A KD= 

H h KD= 

X — X KD= 

❖ e>  KD= 

4 ^ KD= 
X — X KD= 


KD= 

KD= 

KD= 


Fig.  30.  Magnitude  of  for  a wire  of  radius 
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Fig.  33.  Magnitude  of  t for  a wire  of  radius  kr 


TgCo.kd.kr) 


The  contour  of  integration  in  ri54)  is  shown  in  Fig. 36  which  assumes  that 
the  medium  is  very  slightly  lossy  so  that  k is  complex.  Branch  cuts  in  the 
plane  are  chosen  along  the  loci  Im  k^  = 0 and  k^  is  assumed  to  be  in  the 
sheet  lmkp<0  so  that  the  radiation  condition  is  satisfied  for  each  component 
cylindrical  wave  in  the  integral  representation  (154).  If  z is  positive,  the 
contour  then  may  be  deformed  around  the  branch  cut  as  shown  in  Fig. 37  (where 
k has  also  been  allowed  to  become  real)  and  the  integral  becomes 


PATH  OF 
INTEGRATION 


BRANCH 

CUTS 


Branch  cuts  and  contour  of  integration  in  the  plane 


DEFORMED / 
PATH 


Deformation  of  the  contour  around  the  branch  cut 


= Ij^C-z),  z > 0 

(157) 

The  first  two  terms  on  the  RHS  of  (156)  and  (157)  represent  the  TENl  modes 
excited  by  the  electric  and  magnetic  dipoles,  respectively,  and  arise  from 
the  branch  point  contribution  from  the  deformed  path  in  Fig.  37  . The  remaining 
two  integrals  in  each  expression  represent  the  current  due  to  the  propagating 
and  the  evanescent  spectrum,  in  that  order.  Noting  that  the  integrands  for  the 
evanescent  fields  are  bounded  by  some  number  C,  say,  we  can  bound  those  inte- 
grals by  the  integral 

Ce‘“^da  = (158) 

0 

The  integrands  for  the  propagating  fields  are  of  bounded  variation  and  hence, 
by  the  Riemann-Lesbegue  lemma  [23]  are  0{l/z)  as  well.  Hence,  for  large  z, 

[23]  E.  T.  IVhittaker,  G.  N.  Watson,  A Course  of  Ntodem  Analysis,  New  York: 
Cambridge  Press,  1973,  p.  172. 
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APPENDIX  A 

LUMPED  ELEMENT  TRAVELING  WAVE  SOURCE 

In  many  network  computations  it  is  convenient  to  use  scat- 
tering parameters,  defined  in  (ref.  Al) , (ref.  A2) 


2/R 


Ok 


(Al) 


b 


k 


(A2) 


In  the  above,  and  Ij^  are  the  actual  (not  normalized)  voltage 
and  current  at  port  k of  an  N-port,  while  = R^j^  + j^Ok 
normalization  impedance  at  port  k.  The  asterisk  denotes  a com- 
plex-conjugate number.  When  scattering  representation  of  net- 
works is  used,  it  is  convenient  to  replace  the  voltage  and  cur- 
rent sources  with  corresponding  traveling-wave  sources.  It  seems 
that  no  convenient  circuit  elements  have  been  found  for  that  pur- 
pose. Penfield  and  Rafuse  (ref. A3)  have  used  a symbol  reminis- 
cent of  a voltage  source  in  combination  with  a directional  coupler 
such  as  in  Figure  Al.  However,  this  symbol  represents  a pure 
traveling-wave  source  only  in  the  limiting  case  when  the  coupling 
tends  to  zero  and  the  voltage  source  tends  to  infinity.  For  any 
finite  value  of  the  coupling,  a fraction  of  the  energy  from  the 
main  guide  will  be  unavoidably  absorbed  in  the  matched  termination 


Al.  D.  C.  Youla,  "On  Scattering  Matrices  Normalized  to  Complex 
Port  Numbers,"  Proc.  I.R.E.  Vol.  49,  p.  1221,  July  1961. 

A2.  M.  T.  Carlin,  A.  B.  Giordano,  Network  Theory,  an  Introduction 
to  Reciprocal  and  Nonreciprocal  Circuits,  Englewood  Cliffs: 
Prentice  Hall,  1964,  pp.  526  and  144. 

A3.  P.  Penfield,  Jr.,  R.  P.  Rafuse,  Varactor  Applications, 
Cambridge:  M.I.T.  Press,  1962,  p"!  26"! 
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A convenient  source  of  the  traveling-wave  to  be  proposed 
here  is  a combination  of  the  voltage  source  and  the  current 
source  such  as  shown  in  Figure  A2(a).  To  verify  this,  write  the 
Kirchhoff  laws  for  this  circuit: 

II  + I2  Is  “ ° 

Vi  - V2  + Vs  = 0 (A4) 

Take  the  normalization  impedance  for  the  port  1 to  be  and 
for  the  port  2 to  be  2^2*  Then  using  (1)  and  (2)  obtain  the 
following: 


A convenient  choice  of  the  normalization  impedance,  is  as 
follows : 

^01  ■ ^02  ■ ^0  ■ ''0  * ^=^0 

When  such  a choice  is  made  the  diagonal  terms  in  (A5)  become  zero 
and  the  two  port  in  Figure  A2(a)  becomes  an  allpass: 


Znl.-V 

Os  s 


(A7) 


Furthermore,  the  source  of  the  outgoing  wave  at  port  1 can  be 
made  equal  to  zero  by  choosing 


(A8) 


Therefore,  under  constraints  (A6)  and  (A8) , the  circuit  in  Figure 
A2(a)  becomes  a true  source  of  an  outgoing  wave  emerging  from 
port  2.  Matrix  equation  (AS)  reduces  to 


where  the  outgoing-wave  source  b 


one  chooses  a real  normalization 


so  that  the  corresponding  outgoing-wave  source  becomes 


From  (A9)  it  is  apparent  that  any  wave  a2  impinging  upon  port  2, 
passes  through  the  circuit  in  Figure  A2(a)  unaffected,  and  emerges 
at  port  1 as  an  outgoing  wave  bj^.  On  the  other  hand,  it  follows 
from  (AlO)  that  wave  a^^,  incident  at  port  1,  passes  unaffected 
through  the  circuit  and  emerges  at  port  2.  In  addition  to  a^^, 
a source  wave  b2g  also  emerges  at  port  2.  This  source  wave  is 
independent  of  b^  or  02*  Hence,  it  has  been  shown  that  a source 
of  a traveling  wave  is  a twoport  consisting  of  a voltage  source 
and  a current  source  such  as  in  Figure  A2(a).  These  two  sources 
must  be  related  by  (A8)  and  the  normalization  impedances  must  be 
selected  according  to  (A6)  fcr  the  complex  normalization,  or 
according  to  (A12)  for  the  real  normalization. 


» 

f 

4 

I c 


To  simplify  the  circuit  diagrams,  the  traveling-wave  source 
between  ports  1 and  2 could  be  represented  by  a single  symbol 
such  as  in  Figure  A2(b).  Note  that  the  symbol  represents  a two 
port  which  is  a three- terminal  device.  The  value  of  the  source 
b2g  is  given  by  (All)  for  complex  normalization  and  by  (A13)  for 
real  normalization. 
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APPENDIX  B 


MATRIX  ALGEBRA  IN  DIRAC'S  NOTATION 


The  presented  material  is  an  abridgement  from  (ref.  Bl) 
and  (ref.  B2).  A column  vector  in  an  N-dimensional  space  is 
denoted  by 


A square  matrix  in  an  N-dimensional  space  is  denoted  by 


Matrix  A is  an  operator.  When  it  operates  on  vector  lx> 
result  is  a new  vector  lv> 


Using  the  rule  for  matrix  multiplication,  one  finds  that  the  k 
component  of  the  vector  |v>  is 


Bl.  B.  Friedman,  Principles  and  Techniques  of  Applied  Mathematics 
New  York:  Wiley,  1956,  pp.  1-33. 


B2.  A.  Messiah,  Quantum  Mechanics,  Vol  I,  Amsterdam:  North 
Holland,  1965,  pp.  162-1^9. 


The  transpose  conjugate  of  column  vector  |x>  is  denoted  by  <x| 
and  it  represents  a row  vector  as  follows 


where  * denotes  a complex-conjugate  value.  The  transpose  con 
jugate  of  a square  matrix  A is  denoted  A 


The  transpose  conjugate  of  a product  of  several  matrices  is 
obtained  by  reversing  the  order  of  multiplication  and  by  taking 
the  transpose  conjugate  of  each  individual  matrix  in  the  product 


The  same  rule  applies  if  one  of  the  matrices  is  a vector  (*N  x 1 
matrix) : 

(A  B |x>)'^  = <x|a'^ 


The  result  of  the  above  operation  is  a row  vector 
which  is  equal  to  its  own  transpose  conjugate: 


A matrix 


is  said  to  be  Hermitian.  Instead  of  distinct  elements,  such 
a matrix  contains  only  ^ N(N+1)  distinct  elements.  The  elements 
located  symmetrically  across  the  main  diagonal  are  complex  conju- 
gates of  each  other.  Also,  the  elements  on  the  main  diagonal  are 


all  real. 

Frequently,  a product  of  a row  vector  and  a column  vector 
is  to  be  computed  as  follows 


<x|y> 


(Xj  X2 


* 

Xn) 


^1 

^2 


N 


* 


3'J 


(B5) 


The  result  of  this  operation  is  a scalar  (=  a complex  number) , 
and  (B5)  is  called  scalar  product.  The  same  two  vectors  |x> 
and  |y>  can  be  used  to  form  a scalar  product  in  the  reverse  order 


<y|x>  - (yj  y* 


I J 


N 


j*l  ^ ^ 


The  result  is  equal  to  the  complex  conjugate  of  CBS).  Thus, 


<x|y>  = <y|x> 


If  one  of  the  vectors  in  a scalar  product  has  been  obtained 
by  the  linear  transformation  such  as  in  (Bl)  above,  the  scalar 
product  takes  the  following  form 

<v|y>  ■ <x|A*|y>  (B6) 

Again,  the  result  of  this  operation  is  a scalar.  The  right-hand 
side  in  (B6)  may  be  interpreted  as  a scalar  product  of  the  row 
vector  <x|A*  with  the  column  vector  |y>.  Alternatively,  (B6) 

98 


may  be  interpreted  as  a scalar  product  of  the  row  vector  <x| 
with  the  column  vector  A*|y>.  Computing  the  product  in  either 
way  gives  identical  results,  because  a matrix  multiplication  is 
associative : 

A (B  g » (A  B)  C . 


When  the  order  of  multiplication  in  (B5)  is  such  that  a 
column  vector  multiplies  a row  vector: 


|y><x| 


'^1  ' 
^2 


* * 

(Xi  X2 


Xn) 


the  result  is  a square  matrix; 


|y><x| 


^1^1  ^1^2 


Vl  V2 


* yA 


(B7) 


This  is  a special  class  of  a square  matrix,  defined  by  only  2 N 
distinct  numbers  (y^  through  yj^  and  x^  through  Xj^) . 

One  of  the  advantages  of  writing  the  matrix  equations  in 
Dirac's  notation  is  that  the  notation  itself  clearly  indicates 
what  type  of  the  result  one  should  expect.  For  example,  the 
expressions  like 

<x|A|x>  , 0<c|x>  , <x|A*B|y> 


represent 

like 


scalars.  Another  class  of  expressions  which  looks 


A|x> 


> 


!x><y|c>  , I a^|x^> 
n-l 

represent  vectors.  Finally,  the  operations  like 

ABC*  or  A]x><b|B|y><c| 

will  each  result  in  a square  matrix.  For  the  sake  of  clarity 
one  should  specify  the  following  was  assumed: 

6,  Ujj  • • • are  scalars, 

l^n^*  |y>^»  |b>,  |c>  • • • are  vectors, 

A,  B,  C • • • are  square  matrices. 
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APPENDIX  C 


MULTICONDUCTOR-LINE  FORMULATION  BY  SIMULTANEOUS 
DIAGONALIZATION  OF  TWO  MATRICES 

Multiconductor  transmission  line  (abbrev.  MTL)  consists  of  N 
conductors  and  a reference,  or  a "ground"  conductor.  The 

voltages  of  each  conductor  with  respect  to  reference  are  arranged  in 
voltage  vector  |V> 

|V>  = 

and  the  currents  are  similarly  arranged  in  current  vector  |l>.  The 
length  coordinate  of  the  MTL  is  z,  while  the  cross-section  coordinates 
are  x and  y.  The  reference  directions  of  the  currents  and  voltages  on 
any  of  the  conductors  are  such  that  real,  positive  values  of  and  I^ 
result  in  a power  flow  in  the  positive  z direction  (see  Fig.  Cl).  On 
the  lossless  MTL,  the  sinusoidal  steady-state  voltages  and  currents  are 
related  through  the  following  equations  (see  e.g.  references  [Cl]  through 
[C8]) 

= -ja,Lll>  , (Cl) 

^|I>  * -ja)KlV>  . (C2) 

L is  the  inductance  matrix,  and  K is  the  electric  induction 
coefficient  matrix.  Often  a s>’mbol  C is  used  to  denote  the  induction 
coefficient  matrix.  Unfortunately,  some  recent  publications  associate 
the  erroneous  term  "capacitance  matrix"  with  the  synibol  C (see[C8]  and 
[C9],  although  the  proper  term  should  be  "electric  induction  coefficent 
matrix"  [CIO] . To  avoid  the  possible  confusion,  symbol  K has  been  used  here. 
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Only  the  lossless  lines  without  magnetic  materials  and  with  iso>* 
tropic  dilectrics  will  be  considered  here.  When  the  dielectric  is 
hoi.iogeneous , the  well-known  TBl  propagation  occurs,  in  vdiich  all  the  waves 
propagate  with  the  velocity  of  light.  Such  ^^^L  are  commonly  treated 
in  terms  of  impedance  parameters  [Cll] , [C8],  image  parameters 
scattering  matrix  parameters  [C13],  and  chain  parameters  [C14] . 

When  the  dielectric  between  conductors  is  not  homogeneous,  the  waves 
on  NTTL  are  HBl  [C15],  i.e,  both  the  electric  and  the  magnetic  vectors 
acquire  a longitudinal  component.  However,  the  dcaninant  HBl  modes  do 
not  have  a cutoff  frequency,  and  at  very  low  frequencies,  their  longi-. 
tudinal  field  conqxjnents  are  negligibly  small.  For  this  reason  these 
dominant  HBl  waves  are  called  "quasi-TEM"  waves,  which  are  the  subject 
of  the  present  report. 

The  formulation  of  voltages  and  currents  on  ^T^L  in  terms  of  voltage 
and  current  eigenvectors,  as  introduced  by  Amemiya  [C6]  and  Marx  [C7], 
is  well  suited  to  the  treatment  of  transients  as  well  as  steady  state 
waveforms.  Their  analysis  starts  frcxn  the  decoupled  equations,  vdiich 
are  obtained  from  (Cl)  and  (C2) : 


- ,!V>  = -oj^LK  |V> 

CC3) 

dz^ 

2-  ,|I>  = KL  |I> 

(C4) 

dz^ 

In  the  present  report,  the  eigenvector  formulation  of  Amemiya 
and  Marx  will  be  extended  by  using  a simultaneous  diagonalization 

of  two  matrices  [C8],  [C16],  [C17].  For  this  purpose,  it  is  convenient  to 
rewrite  (C3)  and  (C4)  as  follows: 
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L'^  ~ o|V>  = -u^K|v>  , 

~ dz^ 

2 

K'^  - 7|I>  = '<i)^L|l>  . 
dz^ 


(C5) 

(^6) 


The  method  of  simultaneous  diagonal izat ion  of  two  real,  symmetric 
matrices  [C16]  presents  a convenient  way  of  solving  these  tvro  equations. 
The  method  can  be  applfed  only  when  one  of  the  two  matrices  is  positive 
definite.  Since  both  L and  matrices  represent  the  stored  energy  of  a 
passive  multiconductor  network,  they  are  both  positive  definite  and  thus 
qualify. 

The  first  step  in  the  method  is  to  find  the  eigenvectors  lx^>  and 
eigenvalues  of  the  matrix  L 

L’^|x.>  = X.|x.>  . (G7) 

The  corresponding  eigenvalue  equation  is 


det(L'^-  XU)  = 0 . 


(C8) 


where  U is  the  identity  matrix.  It  is  convenient  to  make  eigenvectors 
orthonormal : 


<Xj^IXj>  = (Kroneker  delta). 


CC  9) 


Next,  a square  matrix  G is  formed,  having  eigenvectors  |xj^>  for 
its  col’nnns 

G=  (|x^> |3^>)  . (CIO) 


G is  orthogonal  real  matrix  which  diagonalizes  L ^ as  follows 

fX,  0 


g'"l'^g 


•X, 


■N 


= ' A 


(Cll) 


i 
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where  superscript  + denotes  a transpose  conjugate  Since  L'^  is 
positive  definite,  are  positive  real.  Therefore,  it  is  meaning 

ful  to  define  a square-root  matrix 


> 


(C12) 


such  that 

A‘**A  A’**  = U , 
with  U being  an  identity  matrix , 

Next,  a new  variable,  |y>,  is  defined  in  terms  of  |V>  : 


|V>  = G A**|y>  . (C13) 

When  this  change  of  variables  is  introduced  in  (C5) , the  following 
equation  is  obtained 

^ Jy>  = B!y>  , (C14) 

dz^ 

where  B is  a symmetric  real  matrix 

B = A'Vk  G A'^  , • (CIS) 

where  (Cll)  has  been  used  to  eliminate  matrix  L. 

Thus,  B can  be  diagonalized  by  an  orthogonal  transformation  matrix, 
which  will  be  here  denoted  by  The  eigenvalues  of  B are  all  real 
positive.  It  comes  out  that  they  have  a dimension  of  inverse  velocity 
squared: 

B|ti>  = ^2  l^i>  CC16) 

where  |t^>  are  eigenvectors  of  B.  Transformation  matrix  consists  of 


orthonormal  eigenvectors  |tj^>  as  its  columns; 


fiv 


(C17) 


2^  diagonalizes  B as  follows 


2*B2  = 


(CIS) 


(Cl 4)  can  be  thus  diagonalized  by  the  following  change  of  variables 


|y>  = £Ll^> 


(C19) 


Using  (C19)  in  (C14)  one  obtains  the  diagonal  equation 


bJ  0 ' 

» 

• 2 

0 


(C20) 


where 

7 2 

% for  i = 1,2,....N  . 

^ vf 

1 

(C20) is  a decoupled  system  of  differential  equations 


(C21) 


— , C.  = -6?e-  for  i = 1,2,....N  , 
dz^  ^ ^ ^ 

which  has  exponential  solutions 


^i*  ^i 


, e ^ + 5,  e ^ . 


(C22) 


(C23) 


These  are  the  normal  modes  (or  fundamental  modes)  of  the  multiconductor 
transmission  line,  each  of  them  propagating  with  its  own  propagation 
constant  B^j^  in  the  positive  or  negative  z direction.  The  total  vector 


I 


|C>  is  then 


N -jB.z  jB.z 

|C>  = e ^ * qe  ^ )|u.>  , 

i=l  1 1 1 


(C24) 


where  |uj^>  denotes  a unit  column  vector  having  all  zero  components  except 


the  i-th: 


lu.>  = 


i-th  row 


(C25) 


Vector  lc>  thus  represents  the  state  of  the  multiconductor  transmission 


line  in  terms  of  individual  modes,  the  component  describes  the  magni- 


tude of  the  i-th  mode  at  some  position  z along  the  line.  Each  mode  has 
in  turn  certain  magnitude  of  voltage  and  current  on  each  of  the  conductors 
of  the  nulticonductor  transmission  line.  The  voltages  on  individual 
conductors  can  be  now  deteimined  by  transforming  the  coordinates  |5> 
back  into  |V>,  where  the  transformation  is  specified  by  (Cl 3)  and  (C19) 


1V>  = G a'^  g |C> 


(C26) 


The  voltages  on  individual  conductors  are  now 


N -jB.z  _ jB-z  , 

|V>  = I (C.  e * ^ e ^ ) GA  ^|t.> 
i=l  ^ 1 1 


(C27) 


In  the  above  (Cl  7)  was  used  to  write  Q^|u^>  = |t^>  , The  current  vector 
may  be  now  obtained  from  |V>  by  using  (Cl): 


N B.  . -jS.z  jB-z 


(C28) 
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The  power  carried  along  the  multiconductor  transmission  line  is 

P - I Re  <V|I>  (C29) 

where  <\^  denotes  a transpose  conjugate  of  |V>,  and  Re  denotes  the  real  part. 
For  simplicity,  consider  only  the  waves  propagating  in  positive  z direction 
• 0 for  all  i) . Then,  (C29)  gives  the  following  simple  expression  for 
the  power  transmitted  in  the  positive  z direction: 


(C30) 


1 ^2 

The  i-th  mode  carries  the  power  j [C  | /v^.  In  analogy  with  the  two-wire 
transmission  lines  (which  e.g.  consist  of  a single  conductor  and  a shield), 
the  forward-wave  scattering  amplitude  of  the  i-th  mode  will  be  denoted 
by  a^,  and  the  reverse-wave  amplitude  with  b^: 


(C31) 


1 2 

so  that  the  forward- traveling  wave  carries  the  power  , and  the 

1 2 

reverse- traveling  wave  the  power  jl^il  • terms  of  scattering 
amplitudes,  CC27)  becomes 


|V> 


N -jB.z  j$. 
y (a.e  ^ +b.e  ^)|<t).> 

i=l  ^ 1 '^1 


z 


(C32) 


vdiere  the  voltage  eigenvector  is  defined  as 


G A'**|ti>  . (C33) 

The  current  vector  1 1>  from  (C28)  may  be  expressed  in  terms  of  the  same 
voltage  eigenvectors  as  follows: 
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N -jB.z  jB.z 
|I>  - I (a.e  - b.e  J i L 1^ 
i»l  i 


(C34) 


Alternately,  one  can  introduce  the  current  eigenvectors  as  follows: 


and  represent  the  current  vector  as 


(C35) 


N 'je.z  iBiZ 
|I>  = I (a.e  - b.e  ) 

i=l  1 1 1 


(C36) 


If  (C36)  is  substituted  in(C23,  an  alternative  expression  for  lV> 


is  obtained 


N -jB.z  jB.z 
|V>  = I (a.e  ^ + b.e  M 7 K 

i=l  ^ ^ '^i  ^ 


Comparison  with  (C32)  yields 


(C37) 


Normalized  voltage  eigenvectors  |(J)^>,  defined  by(C33),  are  also 
eigenvectors  of  the  matrix 


LK|(J>.>  = i 2 » 

J*  \T  ^ 1 


(C38) 


what  can  be  verified  by  using  (C33)  and  (C16) . Similarly,  is  the 
eigenvector  of  the  matrix  IC  L 


(C39) 


Voltage  and  current  eigenvectors  form  the  bi -orthonormal  set 


(C40) 


Frcxn  (C37)  and  (Cll)  the  explicit  expression  for  current  eigenvector 


is  obtained 


Following  the  notation  frc3in[C6]  and[C7],  define  My  as  a matrix  which 
contains  the  voltage  eigenvectors  as  its  columns: 


Similarly  M,  is  a matrix  containing  all  the  current  eigenvectors 


An  arbitrary  forward  traveling  wave  is,  from  CC32) 


Take  the  scalar  product  from  the  left  with  <ij; 


Thus,  the  expansion  (C44)  can  be  also  written  as 


Since  the  bracket  is  a scalar,  one  can  place  it  behind  |(J).> 


The  above  is  true  for  arbitrary  |V^>,  thus  the  following  identity  is 


obtained 


which  can  also  be  written  as 


Similarly,  one  obtains 


Ml  My  “ U 


M V Mj  = U . 


(C46) 


(C47) 


When  the  propagation  is  assumed  to  be  in  the  forward  z direction,  as  in 
(C44),  the  current  vector  for  an  arbitrary  forward- traveling  mixture  of 


modes  is 


ll£>  * I a.e’j®^^l,j;.>  . 

i=l  ^ 1 


(C48) 


By  multiplying  from  the  left  with  < ())j  | one  obtains 

-jB.z 

a^e  « <<^llj>  . 

Conqjarison  of  this  equation  and  (C45)  yields 

N 

|V£>  * . 

Again,  the  bracket  is  a scalar,  thus  it  may  postmultiply  : 

N 

|Vf>  = n<l>i>  <<{>illf>  * ^|l£> 

The  matrix  relating  |V£>  and  ll£>  is  the  characteristic  iii5)edance  matrix 


(C49) 


This  can  be  written  as 


lo  ' !lv  tC  ■ ’ 


(C50) 


vihere  also  (C46)has  been  used.  Similarly  one  finds  the  characteristic 


admittance  matrix  as 


Multiplying  (C49)  fran  the  right  by  and  using  the  biorthogonality 
property  (C40)  one  obtains 


which  has  the  inverse 


Characteristic  impedance  and  admittance  matrices  can  be  also  related 
through  matrices  L and  K.  The  relationship  can  be  obtained  by  sub- 
stituting (C35)  into  CC42)  and  (C43) 


where  the  diagonal  velocity  matrix  is  defined  as 


When  (L37)  is  substituted  into  (C42)  and(C43)  it  follows 


If  these  equations  are  now  substituted  in  (C51)  the  following  relation  is 


obtained 


For  each  mode  i,  there  is  a fixed  ratio  of  the  voltage  and  current 
on  conductor  j . This  ratio  is  also  called  a characteristic  impedance  of 
the  conductor  j for  the  mode  i,  .denoted  zl.  [C18].  These  characteristic 


impedances 


(scalars)  are  not  to  be  confused  with  the  elements  of  the  characteristic 


admittance  matrix  7^,  If  the  coordinates  of  the  eigenvectors  are  denoted 


by  superscripts,  such  as 

\^i> 


then,  for  the  mode  i,  the  characteristic  impedance  of  the  j*th  conductor  is 


' <UjlV 


(C58) 
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